When does coarsening occur in the dynamics of one-dimensional fronts ? 
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Dynamics of a one-dimensional growing front with an unstable straight profile are analyzed. We 
argue that a coarsening process occurs if and only if the period A of the steady state solution is 
an increasing function of its amplitude A. This statement is rigorously proved for two important 
classes of conserved and nonconserved models by investigating the phase diffusion equation of the 
steady pattern. We further provide clear numerical evidences for the growth equation of a stepped 
crystal surface. 
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Introduction - One of the most fascinating feature of 
nonequilibrium systems is their ability to build up a spa- 
tial (generally) ordered pattern from an initially structur- 
less state when some control parameter reaches a critical 
value. Typical examples can be found in hydrodynam- 
ics, crystal growth, the Turing reaction-diffusion systems, 
and so on These pattern forming systems can be 
broadly classified into two important categories: (i) the 
first category includes the systems which select a typical 
pattern length scale A which is fixed for a given value of 
the control parameters. In that case the whole pattern 
needs not to be steady in time but can undergo oscilla- 
tions, or even exhibit chaotic dynamics (a prominent ex- 
ample is the Kuramoto-Sivashinsky dynamics Q). What 
matters is the selection of a length scale, (ii) The sec- 
ond category corresponds to the situation where there is 
a perpetual coarsening, in that the wavelength A(i) in- 
creases without bound in the course of time. Interrupted 
coarsening may also happen for special situations (see 
below). One must exclude the trivial case where the bi- 
furcation occurs at a wavenumber q = q c and the band 
A.q of active modes extends from q c — 8i to q c + 5 2 - If 
Aq = S1+S2 is not too large, the system selects, of course, 
a pattern with a given intrinsic lengthscale. Examples 
include Turing patterns, Rayleigh-Benard convection..., 
close enough to the instability threshold. Our discus- 
sion focuses on situations where there is no lower bound- 
ary on the active wavenumbers (arbitrarily long wave- 
lengths are active) . It is only at very large distance from 
the instability threshold that, for example, the Turing 
and Rayleigh-Benard systems exhibit nontrivial dynam- 
ics (chaos, coarsening of turbulent plumes...), and this is 
precisely due to the fact that higher and higher harmon- 
ics become active. 

A central question in nonequilibrium pattern forming 
systems is to identify general criteria (if any) for the pres- 
ence of coarsening or not without resorting to a system- 



atic time-dependent calculation. We put forward a cri- 
terion for some classes of nonlinear equations which is 
uniquely based on the behavior of the steady-state solu- 
tions. This criterion is derived on the basis of the anal- 
ysis of the phase diffusion equation drip — Ddxx4>, 
tp being the phase. A negative sign of the phase diffu- 
sion coefficient, D, is therefore a signature for the sys- 
tems that undergo coarsening 3J. The second and more 
important step will be to relate the sign of D with the 
derivative, dX/dA, of the period A of the stationary solu- 
tion with respect to its amplitude A. We will show that 
the sign of D is opposite to that of dX/dA. A coarsen- 
ing process can therefore take place, D < 0, if and only 
if dX/dA is positive. This relation will be proved ana- 
lytically for two classes of nonconserved and conserved 
equations (Eqs. and (fTUf> respectively) and checked 
numerically for an equation which can display different 
dynamical scenarios with varying some physical param- 
eters (Eq. (|14|l ). Finally, considerations based on the 
Lyapunov functional will complement our analysis. 

Nonconserved equations - We begin with the following 
class of nonconserved equations 0]: 

d t u = d 2 u-F(u) = L[u] , (I) 

where dt, d x indicate the time and the spatial derivatives, 
respectively. F(u) is a generic function, playing the role 
of a force. We require that F(u) ~ — u at small u, so the 
trivial solution u(x, t) = is linearly unstable against 
fluctuations of wavelength A > Ao = 2ir. The Ginzburg- 
Landau (GL) equation is recovered for F(u) = —u+u 3 , so 
Eq. will be referred to as the generalized GL equation. 

We can now remark that stationary solutions, uq(x), 
satisfy the condition 

8 2 x u - F(«o) = , (2) 

and uo(x) corresponds to the motion of a particle in the 
potential V(u) = — J duF(u). For a linearly unstable 
profile, V(u) w u 2 jl at small u. 
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The steady-state uq(x) is a periodic solution with the 
spatial periodicity A, uq(x + A) = uq(x). It is useful to 
consider the phase of the pattern, </>, and we shall use 
a nonlinear WKB-like analysis Q] for the derivation of 
the phase equation. For a steady-state periodic solution, 
4> = qx where q = 27r/A. In general the pattern will 
evolve in time, and thus the phase 4>{x, t) is a function 
of space and time. We consider long wavelength mod- 
ulations (which are the most dangerous ones owing to 
translational invariance of the initial pattern), and let e 
denote the smallness of the phase modulation. We in- 
troduce a slow phase ip which is a function of a slow 
spatial variable X = ex. Because we expect any devi- 
ation from the periodic pattern to evolve diffusively we 
anticipate that the slow time T scales as T — e 2 t. The 
slow phase ip(X,T) is related to the fast phase <fi(x,t) 
by <f) — if)/e so that the local wavenumber is given by 
q(X,T) = d<j)/dx = dt/j/dX. In a multiscale spirit we 
have to make the substitutions 



d x -> qd^ + ed x , d t -> e 2 d T (/)d c f, + e 2 d T 



(3) 



We then expand u = uq + eu\ + .... in powers of e. To 
leading order e° we find that uq obeys L[uq] — 0. For the 
generalized GL equation we have 



- F(u ) = 



(4) 



The next order reads L[u\] = f(uo), where L is the 
(Frechet) derivative of L and / is the inhomogeneous 
term arising from phase modulation. For the generalized 
GL equation, L = q 2 d^ — F'(uq), and 

f(u Q ) = dripd^uo - [d^uo + 2qd cf>q u }dxxip ■ (5) 

It is easily shown that, because of translational invari- 
ance, the linear operator L has the nontrivial solution 
d<pU(). It follows that u\ exists only if the inhomogeneous 
term f(uo) is orthogonal to d<piio 6]. This condition re- 
sults in the sought after diffusion equation 



d T ip = Ddxx^ , 
where D, the phase diffusion coefficient, is given by 
d q (q(d<t,u ) 2 ) _ Di 



D 



Do 



(0) 



(7) 



Here above, (...) = (27r) _1 J* ...dcj) is the inner product, 
the denominator D 2 is always positive and the sign of D 
is fixed by the numerator D x = d q (q(d < f,u a ) 2 ) . 

Eq. (@J, which is equivalent to Eq. J2J), corresponds 
to the equation of motion of a particle of unitary mass, 
subject to the force F(uq): uq is the spatial coordinate 
and t = (jy/q is the time. So, we have: 



-M^o) 2 -^, (8) 



where we have introduced the action variable J. 

Finally, if we remember that the period of the os- 
cillatory motion is given by the derivative of the action 
J with respect to the energy (A = dJ/dE), we find: 



27 ~dq~ 



>?_ 

~4ir 2 



dX 
~dE 



(9) 



so that the sign of D\ is opposite to that of dX/dE: 
there is coarsening (D < 0) if and only if the wavelength 
is an increasing function of the energy of the particle or, 
equivalently, of the amplitude A, which is related to the 
energy by E = V(A). 

Conserved equations - We now apply the same pro- 
cedure to a class of conserved equations, which — in 
the same spirit — will be called generalized Cahn-Hilliard 
(CH) equations: 



d t u = -dl[dlu-F{u)] = -d 2 x L[u] 



(10) 



The linear stability analysis of the trivial solution 
u = gives the spectrum uj = q 2 — q 4 and allows to de- 
fine the fastest growing mode, A c = \/2X . We skip the 
calculations to attain the phase equation and just give 
the final expression for the phase diffusion coefficient, 



D = 



q 2 d q (q(d^ ) 2 ) _ q 2 D 1 



(11) 



The q 2 factor at the numerator is distinctly due to the 
conservation law, i.e. to the operator —d 2 in front of L 
in Eq. lfTU|) . The sign of D is fixed by the quantity D\ = 
d q (q(dj,uo) 2 ) as for the nonconserved case. Our criterion 
is therefore proved for the generalized CH equation as 
well. 

Similarly to J — § dT(d T u ) 2 we can define / = 
§ u 2 (t) , where § dr means — in the mechanical analogy — 
the time integral on an oscillation period. Using these 
quantities, D 2 = ((d^,u Q ) 2 ) — XJ/Att 2 and D 2 = (u^) = 
I/X. Finally, the diffusion coefficient D reads 



D = - 



X 2 



J{d E X) 



D = - 



I(d E X) 



nonconserved (GL) models (12) 



conserved (CH) models (13) 



The expressions clearly show that D is negative if and 
only if 8eX > 0. 

The different scenarios - Here above we have estab- 
lished for two important classes of evolution equations 
that the phase diffusion equation is stable (no coarsen- 
ing) or unstable (coarsening), according to the sign of 
dX/dA. Therefore, the following scenarios can be ad- 
vanced (see Fig. 1). (i) If X(A) is a decreasing function 
the system typically develops a profile whose characteris- 
tic wavelength keeps constant in the course of time. This 
happens, e.g., if F = —u — u 3 , i.e. V{u) = u 2 /2 + u 4 /4. 
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(ii) 




(iv) 



FIG. 1: Schematic wavelength-amplitude diagrams for the 
stationary solutions of period A and amplitude A. The four 
scenarios refer to different behaviors for A' (A), as explained 
in the main text. The linear stability analysis of the trivial 
solution u = is represented by the line A — 0: stability (full 
line) for A < Ao = 2tt and instability (dashed line) for A > Ao. 
For the conserved models the flat interface starts to develop 
a profile of period A c = v2Ao (see the big dotted arrows 
close to the A-axis). In cases (ii)-(iv) a coarsening process 
takes place, but in case (iii) it stops when X(A) attains the 
maximum (interrupted coarsening) and the amplitude starts 
growing (see the big dotted arrow close to the maximum). 



(ii) If A (A) is an increasing function the system exhibits 
coarsening. This is what happens in the standard GL and 
CH models: F = —u + u 3 . (iii) If X(A) starts to increase 
and afterwards it decreases the system is expected to dis- 
play coarsening that stops at the maximum of X(A). This 
behavior is midway between (i) and (ii) and it is called in- 
terrupted coarsening, (iv) If X(A) starts to decrease and 
afterwards it increases the system essentially behaves as 
in (ii). 

It is now interesting to consider a growth equation of 
physical interest, which does not fall into the classes we 
have discussed (see Ref. H for a general reference on 
crystal growth dynamics). The equation has the general 
form d t z = -d x j, 

d t z{x,t) = -d x {B{m)+G{m)d x [C{m)d x m}} , (14) 

where m = d x z is the local slope of the front and plays the 
role of u in Eq. i|l(J|) . It is worthnoting that Eq. I|14|l re- 
duces to Eq. ljTU|) if C(m) and G(m) are constant. Equa- 
tions of this kind are often met in nature, such as in 
epitaxial growth of high symmetry as well as vicinal sur- 
faces [H 0|- They are also encountered in the study of 
sand ripples 01 > m problems of spinodal decomposition, 
in the study of many linear defects, and so on. 

It is possible to define the new variable M(m) = 
J dsC(s), so that the stationarity condition j = 



can be rewritten as B(M) + G(M)d xx M = 0. Steady 
states therefore correspond to the solution of the New- 
ton's equation for a fictitious particle, moving in the po- 
tential 



V(M) = / dM 



B(M) 
G(M) 



(15) 



Here we are interested in t he g rowth dynamics of 
stepped, vicinal crystal surfaces [llj by vapor phase de- 
position techniques. With reference to Eq. Ijl4|l . m is the 
local slope of the step and |l2^ : 



B(m) = 
C{m) = 



G(m) = 



1 + /VIT" 



l + m 2 ' ~ v ""' (1 + /3)(1 
l + c(l + m 2 )(l + 2m 2 ) 
(l + c)(l + m 2 ) 3 / 2 ' 



,(16) 



(17) 



where [3 and c are positive, adimensional parameters. 
The former, /3, is the relative strength between line dif- 
fusion and terrace diffusion, and the latter, c, is the 
strength of the elastic coupling between steps. The lin- 
ear spectrum reads ui = q 2 — q A . We now report the 
results of our analysis of A (^4) for the different values of 
the parameters (3 and c (see also Fig. 1). 

If elastic interactions are absent, c = 0, A is always a 
decreasing function of A (regime (i)). If elastic interac- 
tions are weak, < c < 1/(1 + 2/3), X(A) decreases at 
small slopes and increases at large slopes (regime (iv)). 
For strong elastic interactions, c > 1/(1 + 2(3), X(A) is 
always an increasing function (regime (ii)). 

From a dynamical point of view, we therefore expect 
no coarsening for c = and perpetual coarsening for 
c ^ 0. The case without elastic interactions has been 
considered in Ref. for f3 = and in Ref. 9] for any 
13: Refs. 0,13 

clearly show the absence of coarsening and 
the development of a constant pattern length scale with 
A = A c . The case c ^ has been treated in Ref. ( see 
Figs. 1,2 therein) and a persistent coarsening is revealed, 
both for c < 1/(1 + 2/3) (regime (iv)) and for c > 1/(1 + 
2/3) (regime (ii)). 

Recently, Eq. (|14f) has been modified in order to take 
into account anisotropy effects 0|. In this case, A(^4) 
may present a maximum, which leads to regime (iii). As 
expected by our criterion, the coarsening is found [15J to 
interrupt at the maximum of A(^4), where the amplitude 
grows unstably without bound. 

An interrupted coarsening process was also found in a 
discrete model describing the growth dynamics of a one- 
dimensional surface, the so-called Zeno model 16] . The 
surface is described as a sequence of steps which move 
following the Burton, Cabrera and Frank dynamics [l7| . 
What is relevant in our context is that the interruption 
of coarsening at A = A* is related to the disappearance of 
stationary configurations for A > A* 0] (Fig. 1, scenario 
(iii)). 
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Lyapunov functional - It is well known |18( that the 
two classes of nonconserved, Eq. (JIJ, and conserved, 
Eq. ijTUI) . models can be derived by the Lyapunov func- 
tional 19] T[u(x,t)] = J dx [\{d x uf - V(u)], through 
the relations d t u — ~{8T/8u) (GL models) and d t u = 
d x (ST/Su) (CH models). It is also easy to check that 
dT/dt < in both cases, i.e. dynamics proceeds so as to 
minimize T . 

It is trivial that dT/dt vanishes for the stationary con- 
figurations of period A, dtT[u\(x)\ = 0. Rather, we are 
interested here to study the dependence of T on the pe- 
riod A of the steady state. The relation T[u\(x)]/£ = 
( J/A — E) is found, where t is the length of the grow- 
ing front, and the derivative is easily evaluated |20| . 
l- l {dT[u x {x))/dX) = -(J/A 2 ) < 0. This result, together 
with the relation dT j dt < which ensures that dynam- 
ics minimizes T, complements the scenarios described in 
Fig. 1: if A increases with A there is coarsening, and if 
A (A) has a maximum at A = A* coarsening stops at A*. 

Perspectives - Our analysis has focused on equations 
where either coarsening takes place, or the wavelength 
is frozen while the amplitude 'diverges' depending on 
whether X(A) is an increasing or decreasing function. 
This does not exhaust all possible scenarios that occur 
in nonequilibrium systems. A prominent example is the 
Kuramoto-Sivashinsky (KS) equation which is known to 
produce spatio-temporal chaos with a dominant length 
scale (corresponding to the fastest growing mode). It 
would therefore be desirable to see how the steady-branch 
A (A) would behave in the presence of KS-dominated dy- 
namics. We have done a preliminary study j^jj of the 
modified CH equation [22| : d t u = — d x [d x u + u — u 3 } + 
vud x u. For v = this equation reduces to the CH one, 
while for large v an appropriate rescaling leads to the 
KS one. Above a critical value v$ the branch X(A) un- 
dergoes a fold singularity where A ( A) exhibits a turning 
point. We have found that in correspondence of that an 
anticoarsening process appears, that is to say X(t) de- 
creases. Our work is in progress and it aims to clarify 
possible relations between the sign of A' (A) and the type 
of dynamics exhibited by the system. 

A more general puzzling question is whether, for ex- 
tended systems enjoying basic symmetries (translation in 
the plane, x — > x+xq and u — > u+uq (or z — » z + zq), and 
parity) coarsening, chaos or 'diverging' amplitudes with 
a frozen periodicity are the only possible scenarios. Is 
there any simple link between symmetries and the kind 
of dynamics (order or disorder) that a system exhibits 
when it does not succumb coarsening? 

A final question concerns the diffusion coefficient D 
when (perpetual) coarsening occurs, D < 0. In this case, 
the typical pattern length scale A increases in time. Is 
it possible to derive the coarsening law, A(t), from the 
knowledge of |D(A)|? It will be an important task to 



clarify all the above questions and to consider the possible 
extension of our arguments to higher dimension as well 
as to nonlocal equations, as they arise for example in 
solidification and viscous fingering. 

We gratefully acknowledge useful discussions with sev- 
eral members of the Florentine Dynamics of Complex 
Systems group. 
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